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The Maximum Entropy Method is applied to dynamical fermion simulations of the (2+l)-dimensional 
Nambu— Jona-Lasinio model. This model is particularly interesting because at T = it has a broken phase 
with a rich spectrum of mesonic bound states and a symmetric phase where there are resonances, and hence 
the simple pole assumption of traditional fitting procedures breaks down. We present results extracted from 
simulations on large lattices for the spectral functions of the elementary fermion, the pion, the sigma, the massive 
pseudoscalar meson and the symmetric phase resonances. 



1. Introduction 

The QCD spectral function p{uj) carries infor- 
mation not just about particle masses but also 
about resonance widths. Recently, the first se- 
rious attempts to reconstruct spectral functions, 
defined via 

/•oo 

G{t)= / K{t,uj)p{Lu)duj, (1) 
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where the kernel is a free particle correla- 
tor with mass w, from lattice data have been 
made. The inverse problem is numerically ill- 
posed, since p{uj) is a continuous function whereas 
the lattice propagator G{t) is only available for a 
discrete and finite set of points. The procedure 
adopted has been to assume an ansatz for p{ijj) 
based on QCD sum rules [Q. A more general 
solution is to apply maximum entropy method 
(MEM) techniques which extract p{ijj) sub- 
ject to some reasonable requirements that it be 
smooth and stable with respect to small varia- 
tions in the input data G{t), but otherwise make 
no a priori assumptions about its form. The 
theoretical basis of MEM is Bayes' theorem in 
probability theory. Details about the technique 
can be found in A sophisticated technique 
of analyzing data such as MEM is also needed 
to study the complexity of the mesonic spec- 
tral functions in thermal QCD below Tc and at 
the range of intermediate temperatures above Tc 
which is dominated by strongly interacting quarks 



in the mesonic channels. 

To refine these new techniques we started 
an investigation of spectroscopy in general, and 
the MEM method in particular, in the (2-1-1)- 
dimensional Nambu— Jona-Lasinio (NJL) model. 
The continuum Euclidean spacetime Lagrangian 
is: 

Nf 

C = *,(^ + mo + cr + i757r)*i + TT^i'y'^ + tt^)- 

The fields and are four-component spinors 
and the index i runs over Nf fermion flavours. 
The model has an interacting continuum limit 
for a critical value of the coupling 1/5^1 which 
has a numerical value w 1.0/a in the large- A^/ 
limit if a lattice regularization is employed [Q. 
At T = and for mo = the model has a U{\) 
chiral symmetry which is broken spontaneously 
for sufficiently strong coupling > and the 
pion field n is the associated Goldstone boson [|j . 
Despite its simplicity the NJL model has an in- 
teresting particle spectrum. The quark propaga- 
tor is gauge invariant and easily measured, en- 
abling the fermion mass to be monitored as chi- 
ral symmetry is restored. The a and tt mesons 
are represented by bosonic auxiliary fields in the 
code, so the correlation functions in these chan- 
nels, including the disconnected diagrams so ex- 
pensive to calculate in QCD, can be measured 
with high statistics relatively cheaply. The Gold- 
stone mechanism in the NJL model is fundamen- 
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tally difFerent from that in QCD. In the NJL 
model the disconnected diagrams are responsi- 
ble for making the pion light ||]. The connected 
contribution corresponds to a pseudoscalar state 
with a mass which is twice the dynamical fermion 
mass in the \a.rge-Nf limit. Another interesting 
aspect of the model is that in the large- TV/ limit 
the correlator of the a meson has a complicated 
spatial dependence due to a strongly interacting 
fermion-antifermion continuum for k> M^ = 2{a) 
i.e. it has branch cuts in the complex k- 
plane in the ranges [iM„, ioo) and {—ioo, —iMg-], 
instead of a discrete series of poles, corresponding 
to a set of bound states in the channel in question, 
which is the case of QCD. In our simulations we 
check whether next-to-leading order effects make 
the a meson stable by separating a bound state 
pole from the two- fermion threshold. The NJL 
model has also an interesting symmetric phase. In 
this phase the large- iV/ meson correlator does not 
interpolate as a stable massive particle state, but 
rather an unstable resonance and the resonance 
width serves as an inverse correlation length [^ji) . 

2. Lattice simulations and MEM analysis 

The model was discretized by using the stag- 
gered fermion formulation and the simulations 
were performed by using the standard Hybrid 
Monte Carlo algorithm. The lattice action and 
details concerning the algorithm can be found in 
We performed simulations with both large 
Nf = 36 and relatively small Nj = 4. This enables 
us first to compare the MEM results with the 
analytical large-TV/ predictions and then to use 
the method to probe next-to-leading order cor- 
rections. The propagators for the auxiliary fields 
TT and a are given by: 

aW=^(^(x,i')vr(x,t' + t)), (2) 

x,t' 

C.{t) = 5][(a(x, t')a(x, t' + t)) - {af]. (3) 

In order to minimize contamination from un- 
wanted excited states the massive PS meson cor- 
relator was measured by using wall sources. 

Figure |l| shows the spectral functions output 
from MEM for the tt, fermion (f) and massive PS 



meson in the broken phase of Nf = 4 on a lattice 
of size 32^ x 48 and coupling /3 = l/g"^ = 0.55. 
We generated approximately 30, 000 configura- 
tions for each parameter set. 
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Figure 1. Broken phase spectral functions, 
peaks from left to right are for the vr (mo = 
0.005, 0.01, 0.02), f and PS (both mo = 0.01). 
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Figure 2. PCAC scaling relation for f3 = 0.55 data 
sets. 

These particles are all expected to be poles. 
The masses extracted from the MEM anal- 
ysis are in very good agreement with stan- 
dard single-exponential fits. (e.g. for the 
TT in iV/ = 4, /3 = 0.55, mo = 0.01: 1- 
exp=0. 168(6), MEM=0. 162(6) and for mo = 
0.02: l-exp=0.236(6), MEM=0.231(6)). In or- 
der to check that the width of the peaks seen in 
Figure |l| are statistical, an extrapolation to an in- 
finite number of configurations will be performed. 
The binding energy, e^, for the massive PS meson 
is calculated from 2M/ — Mps- For (3 = 0.55 
mo = 0.01: 66 = 0.027 for TV/ = 4 and eb = 0.006 
for Nf = 36. The result is in agreement with the 
expectation that the meson is more tightly bound 
for smaller Nf . Figure || shows that the pion mass 
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obeys the PCAC scaling relationship since the lin- 
ear fit intercepts the y-axis at -0.001(6). 

For the a meson we generated very high statis- 
tics (~ 0(10^)) on a 24^ lattice using the sim- 
pler .Z2-symmetric (rather than J7(l)-symmetric) 
model. It is not clear yet whether the widths of 
the spectral functions shown in Figure || are due 
to the noise in our data set or whether they have 
a physical significance. However, the positions of 
the peaks imply that the sigma is a bound state, 
because e.g., for (3 — 0.70, ujpeak = 0.475(16) is 
significantly smaller than 2My = 0.534(1). This is 
another next-to-leading order effect that we de- 
tected in our MEM analysis. 
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Figure 3. Sigma meson in Z2 model, TV/ — 4. 

In Figure |^ the results obtained for the sym- 
metric phase of the Z2 model {Nf = 4 on a 
32^^ X 48 lattice) are presented. The critical cou- 
pling is (3c = 0.84(1) 1^ and the simulations were 
performed with (3 = 0.92, 1.0, 1.25. The expected 
features (in the large- A^/ limit) are a broad reso- 
nance whose central position and width increase 
with j3. It is not clear yet whether the peaks near 
the origin have a physical significance or whether 
they are finite Lt effects. The correlators in the 
symmetric phase decay very slowly and remain 
significantly non-zero at large t. This behaviour 
is consistent with the power-law form prediction 
of the large- A/ calculation We are currently 
performing new simulations on lattices with both 
smaller and bigger Lt in order to gain a better 
understanding of p{uj) near the origin. 

3. Summary and Outlook 

The preliminary results presented in this report 
show that MEM is a reliable method to under- 
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Figure 4. Symmetric phase scalar for Nf = 4. 



stand the details of the particle spectra of NJL 
and other models. We are currently increasing 
our statistics and also intend to calculate errors 
on the spectral features in the MEM results. Fur- 
thermore, we are also planning to extend our in- 
vestigation to understand the spectral properties 
of the mesons at high temperature and density. 
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